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Lattice QCD at maximal twist 

R. Frezzotti'' and G.C. Rossi^ 

^Dipartimcnto di Fisica, Universita di Roma Tor Vergata and INFN Sezione di Tor Vergata 
Via dcUa Ricerca Scientifica - 00133 Roma (ITALY) 

In this review we discuss the general features of maximally twisted lattice QCD. In particular, we illustrate 
how automatic 0(a) improvement can be achieved and how it is possible to set up a lattice regularization scheme 
where the problem of wrong chirality mixing (be it finite or infinite) affecting the computation of the matrix 
elements of the CP-conserving effective weak Hamiltonian is neatly avoided, while having at the same time a 
positive determinant even for non-degenerate quark pairs. The question of reducing the large cutoff effects that 
appear when the quark mass tends to zero as a consequence of parity and iso-spin breaking in the action is also 
addressed. It is shown that such dangerous lattice artifacts are strongly suppressed if the clover term is added to 
the action or, alternatively, the critical mass is chosen so as to enforce the restoration of parity. 
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1. Introduction and content 

Waiting for new generation computers that 
might allow full-fledged simulations with exactly 
chirally invariant fermions, i.e. fermions obeying 
the Ginsparg- Wilson condition , a viable alter- 
native could be to use maximally twisted Wilson 
fermions |2l3l4l5l5| . possibly coupled with a judi- 
cious choice of the gauge action [7]. 

In this review we wish to outline the structure 
and the properties of maximally twisted lattice 
QCD (Mtm-LQCD) as developed in rcfs. |3l4l5lfij . 
Soon after introducing the idea that to avoid ex- 
ceptional configurations in Wilson fcrmion sim- 
ulations one should introduce quarks in flavour 
pairs and have the Wilson term rotated with re- 
spect to the quark mass term by an axial rotation 
in iso-spin space, it was realized that an especially 
useful choice for that angle is to set it at its max- 
imal value, \uj\ = 7r/2, because in this situation 
0(a) (actually 0(a^'^+^), fc > 0) improvement of 
physical quantities is automatic (sect.|2). 

It was then shown in ^4 that the nice improve- 
ment properties enjoyed by Mtm-LQCD, derived 
for pairs of mass degenerate quarks in [J, can 
be immediately extended to the more interest- 
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ing non-degenerate case, without loosing the pos- 
itivity of the corresponding fermion determinant 
(sect. Ol. This last property is obviously crucial 
if one wants to be able to set up workable Monte 
Carlo-like simulation algorithms for QCD. 

With these ingredients and exploiting the flexi- 
bility offered by the possibility of regularizing dif- 
ferent valence flavours with different value of the 
Wilson parameter, r, it is possible to construct |3] 
a hybrid theory, where sea quarks arc introduced 
as pairs of non-degenerate particles and valence 
quarks are regularized as Osterwalder-Seiler j2| 
fermions, such that no "wrong chirality" mix- 
ing '9 affects the computation of the matrix el- 
ements of the CT^-conserving AS* = 1,2 effective 
weak Hamiltonian (sect.^. 

Despite all these nice features, there remain im- 
portant cutoff effects, originating from the break- 
ing of parity and iso-spin induced by the presence 
of the twisted Wilson term in the action, which 
tend to become large as the quark mass is lowered. 
These lattice artifacts have been discussed both in 
chiral perturbation theory (xPT) JOI-^^li sis well 
as in the language of the Symanzik expansion |S] 
with the conclusion that they can be substantially 
reduced if the clover term "14" is introduced in the 
action with its non-perturbatively fixed csw co- 
eflticient ^Hl or, alternatively, if the critical mass 
is chosen in some "optimal way" (sect.EJ. 
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For lack of space we will not discuss how and 
to what extent the strategy outlined above for 
improvement can be extended to the Schrodinger 
functional formalism |16I17| . Nor we will address 
the very important issue of meta-stabilities p,8l 
I19I20| . that are seen to affect unquenched sim- 
ulations [25 at coarse lattice spacing, and the 
closely related flavour breaking effects visible in 
the value of the square mass difference between 
charged and neutral pions These questions 

and the the present status of quenched and un- 
quenched twisted simulations have been recently 
reviewed in ref. 12'A\ . In the concluding section 
(sect.EJ we will only offer some hints on possible 
ways to overcome (even on coarse lattice spac- 
ings) difficulties associated with the existence of 
metastable phases and flavour breaking effects. 

2. Improvement 

The twisted fermion action for a pair of degen- 
erate quarks in the "physical basis" (where the 
parameter which gives a non-vanishing mass to 
pions, TUq, is real) has the expression 



-mg]^p{x) . (I) 



In eq. we have introduced the following defini- 
tions: 7V = ^7;x(V* + V^) with and V* the 
forward and backward lattice covariant deriva- 
tives, respectively, and Wcr — — J'f V* V^-I-Mcr(r) 
with Mcr{r) the critical mass. We observe that to 
match the r-parity property of the Wilson term 
Mcr(r) is necessarily an odd function of r. The 
form of the fermionic action in the "twisted basis" 
is obtained by performing on the fermion fields 
the axial rotation 



?A , ip ^ ipe 



ji^75 /2 



(2) 



The first observation about Sp (and the one that 
originally led to its introduction) is that the deter- 
minant of the associated Dirac- Wilson operator 
is positive definite for any lo ^ and ruq ^ as 
its expression is seen to be given by 

Vp = det[(£>w + "^9 cos uj)\D^ + niq cos uj) + 



-m^ sin^ UJ 



(3) 



The lattice action is invariant under the (spu- 

rionic) transformations 

1) TZ^ X (r — r) x (rrig — > —mq), where 



7^5 



1p -0' = -"075 



ni 



2) Vci X (r ^ — r) X (rriq —mq), with 



V. 



d ■ 



Uli-x - afi)) 



(0(a;),V(x))-^ e3*^/2 {^(-x),^{-x)) 



(4) 



(5) 



These invariances imply for the lattice vacuum 
expectation values (v.e.v.'s) of (multi-local) oper- 
ators, 0{x) = 0{xi,X2, ■ ■ ■ ,Xn), Xi ^ X2 ^ ■ ■ ■ ^ 

Xn, the relations 



1) (0(x)) 

2) (Oix)) 



(r.nig) 



^ (-1)^'^5[0](0(^)) 

= {-iY^^\o{-x)) 



(-r.-rrig) 



(6) 
, (7) 



where Pfi^^Q] is the T^-s-parity of O and d[0] its 
mass dimension. 

We now prove that i) for generic values of the 
twisting angle, cu, averages (Wilson average - WA) 
of v.e.v.'s computed with opposite values of r are 
free of □(a^'^+i), fc > 0, cutoff effects; ii) for the 
special value \uj\ = 7r/2 averaging is unnecessary. 

These results follow from using the Symanzik 
idea that lattice artifacts can be described in 
terms of the continuum correlators of an effective 
low energy theory renormalized at the scale a~^, 
as well as the symmetry properties of the lattice 
theory and the corresponding ones enjoyed by the 
associated continuum theory. 

i) I a; I 7^ 7r/2 - We write the Symanzik expan- 
sion of the lattice v.e.v. of a multi-local, multi- 
plicatively renormalizable (m.r.) operator in the 
schematic form 



{0{x)) 



aJ2riSA^,r){Oe{x)) 



CS{L0,r){O{x)) + 



(nig) 



+ 0(a2), 



(8) 



where the operators Oi have mass dimension of 
one unit larger than that of O, d[Oe] — d[0]+l. We 
explicitly remark that factors of mq may appear 
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in Og. The lattice invariance Q allow us to write 



{0{x)) 



(-l)'^[oi[cg(^,-r-)(0(-x-)) + 



+0^77° (c.,-r)(0,(-x))' 



+ 0(a2). (9) 



We can now use the invariance "Dd x (m^ ^Tn-q) 
of the continuum theory to bring continuum cor- 
relators computed in —x and — to correlators 
computed in x and iriq. One thus gets 



(Oix)) 



(r,m,) 



Cg(c.,-r)(0(a;)) + 



.^,7a(c^,-r)(0,(:r)) + 0(a2) . 



(10) 



Comparison with JS)) gives the relations 
which imply the WA improvement formula 



{0{x)) 



(UJ) 

(r.niq) 



{0{x)) 



= 2Cg(a;, r)(0(x)) +0(^2) 



(11) 



Exploiting the 7?,5-parity transformation proper- 
ties of O, the previous formula can be rewritten 
in the form of a mass average (MA) 



{0{x)) 



(r\m,) 



+ (_l)^'-«5[Ol(0(a;)) 



= 2Cg(a;, r)(0(x)) +0{a^) 

(m<j) 



(12) 



Note that in the presence of spontaneous chiral 
symmetry breaking the two terms in the l.h.s. 
of p2|l will be in general different even in the 
limit rUq — > 

One can easily convince oneself that the whole 
argument developed above is actually valid for 
any odd powers of a in the Symanzik expansion. 
As a result all 0{a'^''~^^), fc > 0, terms get cancel 



^In the absence of spontaneous chiral symmetry breaking, 
Uke in 2-dimensions, lattice data show improvement at any 
u! even without WA 1251 . 



by taking either the Wilson or the mass average 
of lattice correlators. 

ii) \uj\ = 7r/2 - At maximal twist the lattice 
action is invariant under the further transforma- 
tion (which does not require changing sign to r) 
V xPrfX [niq — !■ —TTiq), whcrc V is ordinary parity 
{xp = (a;o, -x)) 

Uo{x) ^ Uoixp) , (13) 
Ukix) ^ Ulixp ^ ak) , 1,2,3 

This fact allows to prove that (for parity-even op- 
erators) 0(a^'^'*'^), fc > 0, improvement is auto- 
matic. The proof (see Appendix of ref. [S]) closely 
follows the previous line of reasoning and does not 
depend on the r-parity of Mcr- One starts with 
eq. (jHJ taken at, say, uj = tt/2. Using the lattice 
symmetry V x x (rrig — > —vriq), it yields 



{0{x)) 



-it/2 



(^_l^PlO]+d[0] 



(-l)^[01+'^[01(0(-a;p)) 
Cg(r)(0(-xp)) + 



7r/2 



+aJ2vSA^){0,i-xp)) 



O(a^) 



(14) 



In eq. (|14|l P[0] is the parity of O. As before, we 
can act with the corresponding continuum sym- 
metry transformation on the r.h.s. of this equa- 
tion to bring back its space-time argument from 
—Xp to X, obtaining 



{0{x)) 



|7r/2 
(r,m 



C8ir){0{x)) + 



(15) 



-aE'?aW(-l)^[°^l+^[°l(0.(-)> , +0(a^). 

Comparing now eq. H14(l with eq. H15(l we get 

_(_l)P[o.]+P[o]^0(,)^^0(,)^ (16) 

which tells us that, if O is a parity-even oper- 
ator {P[0] = mod(2)), aU 0{a^''+'^) (fc > 0) 
terms in the Symanzik expansion of {0{x))\J^'^ -j 
are multiplied by continuum matrix elements of 
parity-odd operators that hence vanish owing to 
the parity invariance of the continuum theory. 
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3. Non-degenerate quarks 

In this section we want to show that the 
nice properties of Mtm-LQCD, i.e. automatic im- 
provement and positivity of the fermion determi- 
nant, estabhshed for the mass degenerate case, 
can be extended to the non-degenerate case. 

It was proved in ref. @] that the action 

Sf,U^,4>,U]= (17) 
= a'* ^ ■ip{x)['yV - ijnTiWcr + mq + T3eq]ip{x) 

X 

describes the (lattice regularized) gauge interac- 
tions of a pair of non-degenerate quarks with 
renormalized masses TOg^"* = Zp^iriq ± Zg^eq, 
where Zp and Zs are the renormalization con- 
stants of the pseudo-scalar and scalar quark den- 
sity, computed in the massless standard Wilson 
theory. These results follow U from the analysis 
of the chiral WTI's associated to the action IjlTfl . 

We now prove i) automatic improvement and 
ii) positivity of the determinant of the Dirac- 
Wilson operator in (fTTji . 

i) Automatic improvement of v.e.v.'s of parity- 
even operators follows from the invariance of H17|) 
under the transformation Vx-Vdy- {rriq —'mq) x 
(eg — > — Cq) by an argument which exactly paral- 
lels the one developed in sect.|5|for the maximally 
twisted case. 

ii) Positivity of 

<lct[Dnd] = det[7V - i75TiWcr + TOg + TgEg] (18) 

is proved by first noticing that det[£)„d] is even in 
•niq and tq as it follows from the chain of equalities 

det[i?„d] = det[75e"^^"^/'i?„<ie*"^="^/"] - 
= det[Qcr + iTirUq + 75^365] = 

Qcr + £975 ifnq 
iniq Qcr - e<;75 

where we have defined Qcr = 75[7V-|-Wcr] — QIj - 
To proceed it is convenient to further write 

det[Aid] = det[Ql + - el]det[l + 2eqB] , (20) 
B = X7VX = -St , X = {Ql + ml - e2)-i/2 

If m^ — >0, the first factor in the r.h.s. of (pn|l is 
positive. Expanding in SqB, we get for the second 

logdet[l + 2eqB] = tr[log(l + 2eqB)] = 



det 



(19) 



00 ,cy Vk 00 (Af^yi 



fe=i 



= -tr[log(l + 4e2i3tB)]>o. 



(21) 



In the second equality we have used evenness 
in eq and the anti-hermiticity of B. The non- 
negativeness of H21(l is evident as the argument 
of the last logarithm is an operator with a norm 
certainly not smaller than unit. 

We conclude with the following important ob- 
servation. We have proved the positivity (rather 
the non- negativity) of the determinant (|19|l of the 
mass non-degenerate Wilson-Dirac operator un- 
der the assumption that the bare parameters ruq 
and £q satisfy the inequality — > 0. For 
the renormalized masses of any sea doublet this 



bound implies > 



Whether or not 



such a limitation can be a real problem for actual 
simulations is to be seen. We note that in the 
worst case, taking rhq'^^ = rhc ^ 1200 MeV and 
TO^~) = TOs ~ 100 MeV, one gets |^ > 0.85. 

4. Weak matrix elements 

The complicated pattern of "wrong chirality" 
mixing affecting the construction of the renormal- 
izable effective weak Hamiltonian operator on the 
lattice when Wilson fermions are employed 9"26j 
has up to now prevented a reliable, full-fledged 
evaluation of the phenomenologically important 
CT'-conserving weak matrix elements relevant 
for the evaluation of Bk and the amplitudes ap- 
pearing m the famous A/ = 1/2 rule \M)VM\ . 

Mtm-LQCD offers a comparatively simple 
framework where the problem of "wrong chiral- 
ity" mixing can be neatly circumvented leading 
to a cheap computational scheme. As an ex- 
tra bonus automatic improvement, which we dis- 
cussed in the previous sections, is guaranteed 

The idea is to make recourse to a hybrid formu- 
lation in which sea quarks are introduced as max- 
imally twisted mass non-degenerate pairs, while 
valence quarks are taken as Osterwalder-Seiler 
(OS) flavour singlet fermions. This formulation 



*See, however, ref. 1271 for some alternative way-outs if 
the Wilson theory is fully 0(a) improved I28I29I . 
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can be made fully local if for each valence 
quark a ghost field is introduced with the purpose 
of canceling the determinant coming form the va- 
lence quark integration. Provided sea and valence 
quarks are given the same renormalized masses, 
the sector of the theory where only correlators 
with no ghost fields are considered is unitary, with 
"partial quenching" effects due a slightly differ- 
ent sea and valence quarks regularization starting 
only from order a? |33| . 

The key point on which the rest of this sec- 
tion is focused is the fact that killing the effects 
of "wrong chirality" mixing in the lattice data 
of weak matrix elements requires to consider a 
specific regularization of the valence sector of the 
theory which depends on the particular matrix 
element one is willing to compute. In particular 
we will see that it may be necessary to replicate 
certain valence flavours assigning appropriately 
tuned signs to the corresponding Wilson terms. 

Notice that in the Mtm-LQCD hybrid formula- 
tion we are considering parity and flavour break- 
ing effects can only arise from the valence sector. 
The flavour symmetry of valence sector can be 
larger than the direct product of N f U(l)-factors, 
depending on the choice of the values of the Wil- 
son parameters. 

4.1. Bk with no mixing 

In the formal continuum QCD theory with 4 
flavours (QCD4) Bk is defined by the formula ^ 

(Xo|0^^fAA(M)|/^o) - ^^MlFlBM (22) 

^vvtAA = (s7mC?)(s7m^) + (s7M75rf)(s7M75rf) 
and can be extracted e.g. from the correlator 

Ckok = {^K{x)0^%tAA{^)^K{y)) , (23) 
$A' = h^s . (24) 

The key observation of ref. |S] is that the same 
information can be extracted from the correlator 

Ck'QK^ (.^K'S^)'^Q^^V=+AA{mKSy)) , (25) 

^K!^ = d'vl5s'y , = dy-f5Sy , (26) 

^Undefined notations as well as further details on the 
renormalization properties of the different regularizations 
we will discuss below can be found in ref. and will not 
be repeated here for lack of space. 



Qvv+liA — (27) 
= (s„7^d„)(s^7pO + {svlt,l5dv){s'^lt,lzd'^) + 

+ (s«7pO(sWM^f ) + (s,;7/x 750(5^^75^1-) , 

where d'^ and s'^ are replicas of d^ and valence 
quarks with mjj = m'^'^ . This statement can be 
proved by noting that, if QCD4 and the theory 
which for short we will call 4s6v (meaning 4 sea 
quarks - Us, ds, Ss, Cs - and 6 valence quarks - Uv, 
u'y dv, d'y, Sv, Cv) are regularized in the same way 
(say a la GW), the two correlators and 
are equal, simply because they give rise to exactly 
the same Wick contractions. This also implies 
that the renormalization constant of Qyy!jl^^ in 
4s6v is equal to that of Oy^^^iA in QCD4 (if the 
same renormalizalization condition is employed). 
Similarly dyj^Sv and d'^j^s'y have the common 
renormalization constant equal to that of dj^s. 

It is not too difficult at this point to prove that, 
if we take sea and valence quarks of the same 
flavour to have the same renormalized masses and 
we regularize valence quarks as OS fermions hav- 



ing Wilson terms with r„ 



-rj , the 



matrix element (j25|l is not contaminated by un- 
wanted "wrong chirality" mixing effects. In fact, 
the peculiar symmetry properties of the regular- 
ized action prevent the mixing of Qyy!^^^ with 
all the other dimension 6 operators with A5' = 2 
that can be constructed with valence quarks. 

The action of the regularized theory 4s6v^ con- 
structed as described before is invariant (among 
others) under the following transformations 

i) Ex(d„,0 X (m^ ^ ) 

ii) Ex5 (s„,s^) X {ml ^ -m^) 

iii) C X [Ex(d^, s^) x (m^ ^ mf,)] x 
X [Ex5«,s;) X (m^' ^ -m{)] 

iv) 7^5 X (M -> -M). 

where ^ (M is the whole set of mass parameter) 



4' 



4' 



ibqt 
ibqt 



qt 
qt 



-qti5 
-qtib 



'Ps ■■ I ql ix)^-f5-foql{xp) q_l{x)^-qf{xp)-fo-f5 

^To be precise these transformations must be appropri- 
ately extended to the ghost fields associated to each va- 
lence quark in order to make the full action invariant 



6 



and C is charge conjugation. One checks that 
the operator Qyy^j^j^ is even under any of the 
above transformation, while the dimension 6 op- 
erators with which it can potentially mix (namely 
nAS=2 nAs=2 nA-s=2 r>^s=2 

^VV-AA^ ^SS-PP' ^TT ' ^SS+PP' 

Qsp±ls and Q^|=2) are aU odd. 



nAS=2 
^VA±AV^ 



4.2. Kaon decay amplitudes 

A long standing puzzle in low energy hadron 
physics, often referred to as "octet enhancement" 
or "A/ = 1/2 rule" , is represented by the surpris- 
ingly large experimental value of the ratio |3(Jj 



■ tttt) = 



Y(K -> 7r7r)|Aj=i/2 
Y(K 7r7r)|A/=3/2 



400 . (28) 



Though the rate of the A/ — 3/2 K tttt weak 
decays can be reasonably well computed within 
our present understanding of field theory (based 
on O.P.E. and renormalization group arguments), 
theoretical estimates of the A/ = 1/2 amplitude 
give much too small values compared to the ex- 
perimental number |34| . 

The lattice approach seems to be the natural 
framework where a first principle evaluation of 
such amplitides can be carried out. There are, 
however, severe difficulties in the process of es- 
tablishing a proper strategy to accomplish this 
task. First of all, in Euclidean metric the pro- 
cedure necessary to extract the kaon decay am- 
plitudes of physical interest is complicated by IR 
subtleties arising from final (two-pion) state inter- 
actions 35 . To attack this problem new interest- 
ing ideas have been recently put forward 36"37'. 
Here we will only deal with the UV difficulties 
related to the construction of the rcnormalized 
effective weak Hamiltonian operator on the lat- 
tice and show how exploiting the flexiblity of the 
Mtm-LQCD formulation the problem of wrong 
chirality mixing can be brought to an amenable 
solution. For an interesting step forward in this 
direction see also |38| . 

4.2.1. Generalities 

In the Standard Model the decay of the K me- 
son into pions is described to leading order in 
the Fermi constant, Gp, by the CT'-conserving, 
AS = 1 effective weak Hamiltonian, which (in 



the chiral limit) reads [HOI 

7i,^/=i = KdV;;%E^^W^w)a^(A^), (29) 

with VudV*g the product of the appropriate ele- 
ments of the CKM matrix ^. The effective opera- 
tor Ti^ff^^ is obtained after having integrated out 
all degrees of freedom above some energy scale, A, 
with A larger than the charm threshold, but still 
well below the W-boson mass, Mw Consistently, 
the operators in eq. H29|l are rcnormalized at 
the scale /i with nic ^ /U <C Mw, while the Wil- 
son coefficients C±{ii/Myf) carry the information 
about the physics between fj, and Mw. The ex- 
plicit expression of the bare operators correspond- 
ing to (5± is (7^ = 7^(1-75)) 



20±=[(s7»(u7^d) ± {sjf;d)iu-f^u)]-[u^c] . 

To make contact with experimental data it is 
enough to consider the decay of the neutral kaon, 
K°, into either tt+tt" or Tr^TrO states |39I4()I41I38| . 
Owing to the parity invariance of the formal con- 
tinuum QCD4 theory, the relevant amplitudes 
A{K° TT+TT-) and AiK'^ tt°tt°) can be writ- 
ten in terms of matrix elements of the rcnormal- 
ized parity odd operators, Oyj^{fi), whose bare 
expression is 

2CyA = [(«7M'")("7M75rf) ± (s7a<'^)(w7m75u) + 
+ (s7^75u)(w7m'^) =•= (s7M75rf)(u7M")]^[" ^ c] . 
4.2.2. K — > TTTT with no mixing 

The idea of the approach developed in ref. |H] 
is to extend the philosophy employed in sect. 14. ll 
for the computation of Bk to the kaon decay am- 
plitudes. In this case the required pattern of va- 
lence quark replicas will be somewhat more com- 
plicated and one needs to take four replicas of 
both Uv and c„ valence flavours. We immediately 
notice that this entails only a small extra compu- 
tational burden (the evaluation of six rather than 
four valence quark propagators) . 

More precisely it can be shown that any infor- 
mation contained in the QCD4 correlator (meson 
charges are understood) 

QCD4 



C. 



{<P^ix)<PAz)o^Amkiy))\ 



|QCD4 



(30) 

^As usual, the top quark contribution, which is down by 
a factor ©(Vt^yj* /y^^jy^g ~ 10""^), is neglected. 
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can be recovered from the correlator 

C^t]°:, = ($.(x)<I>.(z)Q±^(0)<I>l,(y))r°V (31) 

The niodel 4sl0w is intended as a theory where, 
besides the two doublets of sea quarks (w, d) and 
(s, c), in the valence sector the Uy and Cy flavours 
are replicated four times. Giving them the names 
u^v^ and c^v^ {k — 0,1,2,3), Qy^ must be taken 
to have the expression ® 

Q^. = Q^^i"l + Q^W-iQ^i^l-iQ^i^ (32) 

Qvf-0^^\:zi:i. (33) 

We remark that the interpolating pion and kaon 
fields in eq. besides the quarks d and s, 

should be constructed with only the u'^l replica of 
the up quark. As in sect. 14. ll one can prove that, 
if the two theories QCD4 and 4sl0ti are regular- 
ized in the same way, the correlators H3()|l and (|31|l 
are equal. 

At this point the theory 4sl0t; is regularized by 
taking the signs of the Wilson parameters of OS 
valence flavours to be in the particular relation 

rd = rs = ^^[o] = -r„[i] = r„[2] = -r„[3i = 

= r^io] = -r^ii] = r^i2] = -r^pi . (34) 

and the renormalized masses of (sea and valence) 
quarks of the same flavour are assigned the same 
value. We will give the name 4sl0u^ to this reg- 
ularization of the 4sl0w theory. 

One can prove that in 4sl0ti'^ the opera- 
tors Qyj^ are multiplicatively renormalizable, i.e. 
i) they mix neither with operators of dimen- 
sion smaller than 6, except for (m^ — m,^)(m<; — 
TOd)s75C? which comes into play with a logarith- 
mic divergent coefficient, ii) nor with operators of 
dimensions 6. 

i) We first prove that the pattern of mixing 
with operator of dimension < 6 is as in the con- 
tinuum, namely that the operators of dimension 3 
and 5 [575^, sj^a-Fd] can only appear multiplied 
by the mass factors (771^ — mfj{ms — irid) and the 
operators [sd, sa-Fd] by (tti^ — 7n\){m1 — m^). 
The various mass factors come from the symme- 
tries enjoyed by the regularized theory 4sl0w'^, 
according to the following scheme. 

®From now on not to have much too clumsy notation we 
will drop the subscript „ on valence quark fields. 



• rric — m„ from Y[k Ex(w['^l, c^'^l) x (m„ rric), 

• mg — rriii from Ex((i, s) x C x (mcj ms), 

• mc + mu from [EX5 («["!, ul-^l) x Ex^{u^^\u^^^)x 
•x{mu -ruu)] x [u -> c], 

• ms — rud from V5 x (M ~M), but only in 
association with the operators sd and sa-Fd. 
ii) Proving that the operators Qy^ do not mix 

with operators of dimension 6 can be done by set- 
ting all quark masses to zero {M — 0). Follow- 
ing "Z , we proceed by showing that the symme- 
tries of the regularized theory forbid the mixing 
of anyone of the four operators Qy]f' appearing 
in eq. H32|l with all the other operator of dimen- 
sion 6 having the same unbroken quantum num- 
bers. Let us start by discussing the case of Qyj^*' 

and Qyj^ ■ We first notice that the SU{4) flavour 
symmetry of the massless theory prohibits mix- 
ing between -I- and — operators. Absence of mix- 
ing with operators having tensor structures other 
than VA+AV is proved by identifying the various 
symmetries that forbid them. The result is 

• V5 forbids VV ± AA, SS ± PP and TT. 

• Ex(d, s)xC forbids SP + PS and TT. 

• Ex(d, s) xCx SU{A)y forbids SP - PS and 
VA - AV. 

The analysis of the mixing properties of Qp'^ 
) immediately brought back to the case of 

Qy^A (Qv^a) discussed before, by the observation 
that the change of variables induced by TZ^ on the 
valence quarks u'^l and c'^l (w'^l and c'^l) changes 
sign to and r^[i] (r„[3] and r^[3] ), making them 
equal to r„[o] and r^[o] (r„[2] and r^[2]). 

4.2.3. K ^ TT vifith no mixing 

By the use of soft-pion theorems one can relate 
K TTTT to K ^ TT amplitudes ^^3.9^ in the chi- 
ral limit with the advantage that only three-point 
correlators need be evaluated and no problems 
with final state interactions occur. 

In this approach the relevant target matrix ele- 
ments are, for instance, (7r+, q|C'yy^^^|i4r+, p). 
They can be extracted from the correlators 
c5S4=($.(x)O±^^^^(0)a>t,(y))|QCD4^ (35) 

[(sTp") ("7m^) ± i^lt^d) (ut^m) + (36) 
+ (57^75-") {uif^lsd) ± (57^75^) (-07^75m)]-[u <^ c] . 
Noting that, if QCD4 and 4sl0w are regularized 
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in the same way, the Wick theorem ensures the 
equahty of all contractions, we conclude that the 
same information is contained in 



are affected by dangerous artifacts of relative or- 



(37) 



As before, there is a suitable regularization of 
4sl0w which makes the operators Qyv+AA multi- 
plicatively renormalizable. Besides taking all the 
renormalized quark masses of the same flavour 
equal, one needs to choose the Wilson parameters 
of the OS valence quarks obeying the relation 

Td^-Ts^ r„[oi = -r„[i] = r„[2i = -r„[3] = 

= r^[o] = -r^[i] = r^i2] = -r^oi . (38) 

We notice that the only difference with respect to 
the choice (|34|) is the sign change in front of rg- 
We will call this regularization AslQv^* . 

The proof of this statement follows from the in- 
variance of expectation values under the change 
of fermionic integration variables induced by the 
transformation s — s- s' = 75s, s ^ s' = — S75 on 
the valence strange quark (the associated ghost 
fields should be simultaneously transformed in 
the appropriate way |S]). Under this change 
of variables the action of 4sl0?;^* goes over to 
the action of the 4slOt;^ theory of sect. 14.2.21 
modulo the change of sign of the valence quark 
mass, nis- At the same time Qyy+AA ~^ QvA 
and (rus + ra^jsd (— "t-s + 'mdjs^^d. We 
are thus exactly in the situation discussed in 
sect. 14.2.21 except that now the operators [sd, 
sa-Fd] will appear multiphed by the mass fac- 
tors {ml — m'^){ms+md) and the operators [sj^d, 
S'y5a-Fd] by (m^ — m'f^){ml — m^). 

This argument also proves that the operators 
^vv+AA in QCD4 and Qyy+^^ in 4sl0w^* (can 
be taken to) have equal renormalization constant 
and the same as Q^^ in 4sl0?;^. 

We conclude by stressing that in all the appli- 
cations considered in sect.^ full 0{a^''~^^), fc > 0, 
improvement is ensured by the symmetry of the 
regularized actions under V x T>d x (M —M). 

5. Chirally enhanced cutoff effects 

Although, as we have seen, in Mtm-LQCD odd 
power discretization effects are absent or easily 
eliminated |2I3I4I5| . it turns out that correlators 



der 



,2fc 



fc > 1, which are enhanced by inverse 



powers of the (squared) pion mass, as the lat- 
ter becomes small ff.. In fact, when analyzed in 
terms of the Symanzik expansion, lattice expecta- 
tion values exhibit, as 0, what we will call 
"infrared (IR) divergent" cutoff effects with a be- 
haviour of the form (2fc > h > 1, k, h integers) 



(O) = (O) 



1 + 



-,2k 



(mj)'' 



(39) 



where we have assumed that the lattice correlator 
admits a non-trivial continuum limit. Powers of 
Aqcd required to match physical dimensions are 
often understood in this section. 

We want to prove that artifacts of the type (|39(l 
are reduced to terms that are at worst of order 
a? {a^ / m1)^~^ , fc > 1, if the action is 0(a) im- 
proved a la Symanzik or, alternatively, if the crit- 
ical mass is chosen in some "optimal" way. 

The idea that a suitable definition of critical 
mass exists which can lead to a smoothing out of 
chirally enhanced lattice artifacts or perhaps be 
of help in getting improvement was already put 
forward in the context of lattice xPT in refs. [TU) 
and (m, respectively. 

An important consequence of this analysis is 
that the strong (order of magnitude) inequality 
niq > uAq^j^, invoked in ref. can be relaxed 
to the weaker relation > o^Aq^j^, before large 
cutoff effects are possibly met while lowering the 
quark mass at fixed a. The works of refs. |1(JI11 |. 
and most recently refs. |12ll3j . all based on lat- 
tice xPT, lead to essentially equivalent conclu- 
sions about cutoff effects in pion quantities in 
the parameter region niq > a^Aq^j-,. They also 
yield interesting predictions on the possible Wil- 
son fermion phase scenarios 18 19 and observ- 
ables, when niq is of order or smaller. 

A nice numerical demonstration of the effec- 
tiveness of Mtm-LQCD in killing 0(a) discretiza- 
tion errors and the ability of the optimal choice of 
the critical mass in diminishing the magnitude of 
lattice artifacts at small quark mass can be found 
in refs. |42l43l23j . As for Mtm-LQCD with clover- 
improved quark action, the promising quenched 
tests presented some years ago in |2] have been 
recently extended in |45| down to pion masses of 
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300 MeV or lower, confirming the absence of large 
cutoff effects. 

5.1. Symanzik analysis of "IR divergent" 
cutoff artifacts 

The study of discretization artifacts affecting 
lattice correlators in Mtm-LQCD can be elegantly 
made in the language of the Symanzik expansion. 
We focus here on Mtm-LQCD with two mass de- 
generate flavours. Its fcrmionic action is given 
by eq. with a; = 7r/2. A full analysis of cut- 
off effects beyond 0(a) is extremely complicated. 
Fortunately it is not necessary, if we limit the 
discussion to the terms that are enhanced as the 
quark mass is decreased. 

• The Symanzik LEE A of Mtm-LQCD - The 
low energy effective action (LEE A), S'sym, of the 
theory can be written in the form 



Ss- 



ym 



y 



(40) 



where £4 = 2p*'^(-^ ' ^) + ' ^ + "^"j) ^ ^® 
the target continuum 2-flavour QCD Lagrangian. 
Based on the symmetries of Mtm-LQCD a num- 
ber of interesting properties enjoyed by S'sym can 
be proved which are summarized below. 

1. Lagrangian densities of even dimension, 
£4+2*;, in eq. H40(l are parity-even, while terms of 
odd dimension, ^5+2*; 7 a-re parity-odd and twisted 
in iso-spin space. Thus the latter have the quan- 
tum numbers of the neutral pion. 

2. The term of order a in eq. (|45|) . £5, is given 
(on-shell) by the linear combination 

^5 — S^^SW ib.SW + h,m'^ ^5,m2 + 5^,e ^5,e , (41) 

1/ — — 



^5,e = AQCD-!/;i75T3-0 , 



(42) 



where S^^sw, ^b.m^ and S^^e are dimensionless 
coefficients, odd in r. The operator ^5_e arises 
from the need to describe order a uncertainties 
entering any non-perturbative determination of 
the critical mass and goes together with i^.sw- 
Both t^^sw and ^s.e could be made to disappear 
from 140(1 by introducing in the Mtm-LQCD ac- 
tion the SW (clover)-term ^1] with the appro- 
priate non-perturbatively determined csw coeffi- 



cient II 5J and at the same time setting the critical 
mass to its correspondingly 0(a) improved value. 

3. Higher order ambiguities in the critical mass, 
which will all contribute to >Codd , are described by 
terms proportional to odd powers of a of the kind 

„2fe+l s: — 

0. 05+2fe,e <-5+2fe,e — 

-a^fe+l^g^^fe.e (AQCD)''+'^n5T3V' , fc > 1 . (43) 

• Describing lattice correlators beyond 0(a) - 
We are interested in the Symanzik description 
of the lattice artifacts affecting the correlator 
(O) 1^ , where O has continuum vacuum quantum 
numbers so as to yield a non trivially vanishing 
result as a — > 0. Schematically we write 



{0{x)) = {[0{x) + Aodd0(2:) + AevcnO(a:;)] 



nig 

C30 oo 
£odd — / ^ " ^5+2fe , even — 



,2fe/ 



t4+2fe 



(44) 



(45) 



fc=0 



k=l 



The operators AoddO (AovcnO) have an expan- 
sion in odd (even) powers of a. They can be 
viewed as the operators necessary for the im- 
provement of the matrix elements of O j28ll5| . 
To ensure automatic improvement we shall as- 
sume that O is parity invariant in which case its 
Symanzik expansion will contain only even pow- 
ers of a. 

• Pion poles and "IR divergent" cutoff effects - 
Although a complete analysis of all the "IR diver- 
gent" cutoff effects is very complicated, the struc- 
ture of the leading ones [h = 2k in eq. (EHl) is 
rather simple, as they only come from continuum 
correlators where 2k factors / d'^yCoddiv) are in- 
serted. More precisely the leading "IR divergent" 
cutoff effects are identified on the basis of the fol- 
lowing result inj. 

In the Symanzik expansion of (0)1^^ at or- 
der a^*^ (fc > 1) there appear terms with a 
2fc-fold pion pole and residues proportional to 
|(f}|/:oddk°(0))|2'=, where {n\ and |7r"(0)) denote 
the vacuum and the one-7r'^ state at zero three- 
momentum, respectively. Putting different fac- 
tors together, each one of these terms can be seen 
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to be schematically of the form (£odd — Ol'j)) 

- / 1 \ 2fc -1 cont 

( — ) (^.)2'^X[O;K(0)}2fc] , (46) 

e. = |(17|/:oddk°(0))|5„°f , (47) 

where we have generically denoted by 
M.\0\ {7r°(0)}2fe] the 2A;-particle matrix elements 
of O, with each external leg being a zero three- 
momentum neutral pion. 

Less "IR divergent" cutoff effects (those with 
h strictly smaller than 2fc in come either 

from terms with some extra J d^yCcvcn{y) inser- 
tions or from contributions of more complicated 
intermediate states other than straight zero three- 
momentum pions or from both. In all these cases 
one gets extra powers, not all "accompanied" 
by corresponding l/(m^)^ factors. 

It is important to stress that the appearance of 
pion poles like the ones in H46(l in no way means 
that the lattice correlators diverge as niq — > 0, 
but only that the Symanzik expansion we have 
employed appears to have a finite radius of con- 
vergence (on this point see the remarks in '13^). 

5.2. Reducing "IR divergent" cutoff effects 

Recalling that £odd — + 0{a^), the pre- 
vious analysis shows that at leading order in a 
the residue of the most severe multiple pion poles 
is proportional to |(ri|^5|7r°(0))p'^. It is an im- 
mediate conclusion then that the leading "IR di- 
vergent" cutoff effects can all be eliminated from 
lattice data if we can either reduce £5 to only the 
^5 ,„2 term in (|41() or set to zero. 

• Improving the Mtm-LQCD action - The obvi- 
ous, field-theoretical way to achieve the first goal 
consists in adding the clover term [14115128] to 
the Mtm-LQCD action. In this case lattice corre- 
lators will admit a Symanzik description in terms 
of a LEEA where the operators Is^sw and i^^e 
are absent, and in £5 only i^^m'^ survives. The 
left-over contributions arising from the insertions 
of i^^m^ yield terms that are at most of order 
{am.g/ml.y'' ~ {amq)"^^, hence negligible in the 
chiral limit. It is instead the next odd operator in 
the expansion H4UI) . 0^(7, which comes into play. 

A combinatoric analysis based on the structure 
of the non-leading "IR divergent" cutoff effects 
reveals that the worst lattice artifacts left behind 



in correlators after the "clover cure" are of the 
kind a^{a^/ml)''-\ k > 1. 

• Optimal choice of the critical mass - The al- 
ternative strategy to kill the leading "IR diver- 
gent" cutoff effects consists in leaving the Mtm- 
LQCD action unimproved, but fixing the critical 
mass through the condition 

hm C.(m,)= hm |(r!|£oddk°(0))|^°:* = 0-(48) 

The meaning of (|48() is simple. It amounts to fix, 
for fc > 0, the order 0^*^+^ contribution in the 
counter-term, Mcr'4'^i'y5T3ip'" , so that its vacuum 
to one-7r°(0) matrix element compensates, in the 
limit rUq 0, the similar matrix element of the 
sum of all the other operators making up £5+2k- 

Concrete procedures designed to implement 
condition H48|l in practice were discussed in 
They are all based on the idea of determining 
the critical mass by requiring the lattice corre- 
lator a^EMi^)PHO))L,_ i^o + 0) to vanish 
in the chiral limit {y^ = V'7o"^V' is the vector 
current with iso-spin index 2 and = ipj^^ip 
the pseudo-scalar density with iso-spin index 1). 
In the continuum this correlator is zero by parity 
for any value of rriq. On the lattice the break- 
ing of parity (and iso-spin) due to the twisting of 
the Wilson term makes it non- vanishing by 0(a) 
discretization artifacts, which have the form of a 
power series expansion in 

The important conclusion of the analysis pre- 
sented in [S] is that it is not necessary (nor possi- 
ble) to really go to ruq ^0. It is enough to have 
the critical mass determined by the vanishing of 
the previous correlator at the smallest available 
mg-value, say to™™, provided it fulfills the or- 
der of magnitude inequality to™'" > a^. Under 
this condition, for any rUq > to™'" we will have 
^Trifnq) = O(aTOj). Substituting this estimate 
into the leading "IR divergent" terms, we see that 
they are reduced to finite 0(a^'^) terms. As for 
the sub- leading terms, a non-trivial diagrammatic 
analysis shows that the worst of them, left behind 
after the "optimal critical mass cure" , are reduced 
to only a^(a^/mj)'^"^, k > 1, effects, just like in 
the case where the clover term is employed. 

The results reached in this section show that 
the continuum extrapolation of lattice data 
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should be smooth at least down to values of the 
quark mass satisfying the order of magnitude in- 
equality niq > a^Aqj^p, if the clover term is added 
to the lattice action, as well as if the critical mass 
is adjusted to its "optimal" value. 

5.3. The lattice GMOR relation 

A workable way to numerically estimate the 
minimal value of the quark mass that can be 
safely simulated at fixed a can be obtained by 
considering the behaviour of the charged pion 
mass as a function of ruq. It turns out, in fact, 
that in Mtm-LQCD there are Ward-Takahashi 
identities (WTFs) which take exactly the form 
they have in the formal continuum theory. From 
them a lattice GMOR relation can be derived. 

To see how this works we recall that in Mtm- 
LQCD the 1-point split axial currents, with 
iso-spin index 6=1,2 are exactly conserved at 
TTiq = )2l3j . This implies the validity of the 
WTFs (t± = {Ti±iT2)/2) 



(5°(0)) 4,0 (49) 



5'0=V^V'^-(50) 
After space-time integration, one gets (to^ ^ 0) 
2m,a4^(F±(x)FT(0))|,^,, = -(^"(O)) |,^,, .(51) 

Although, as the WTI ifsT)) itself shows, there is 
no mixing between S'^ and the identity operator 
with a cubically (or a linearly) divergent coeffi- 
cient, there is still room for a quadratically di- 
vergent term proportional to mq. Indeed, the 
l.h.s. of (|51|l can be written as the sum of a piece 
where intermediate states are inserted plus a di- 
vergent contribution of the kind mqja} coming 
from the (integrated) short-distance singularity 
of (P±(x)PT(o)) at X = 0. This term should 
be brought to the r.h.s., thus leading to the sub- 
tracted expression of the chiral condensate. 

We can now repeat on the lattice the argu- 
ment that in the continuum leads to the classical 
GMOR relation. We will work under the assump- 
tion that spontaneous chiral symmetry breaking 
occurs in the limiting continuum theory, i.e. un- 
der the assumption that 

S^-(l]|5"|17)|-"io^0. (52) 



After inserting a complete set of states in the l.h.s. 
of eq. if^ . we get 



(53) 



where we have explicitely written only the contri- 
bution coming from the pion pole. Dots are terms 
due to the intermediate states that stay massive 
as mq 0, as well as terms vanishing with rUq 
faster than linearly. The subscript "sub" is to re- 
mind us that it is the properly subtracted chiral 
condensate that enters this equation. Notice that, 
as expected, the the r.h.s. of eq. (|^ is a finite, 
renormalized, quantity in the limit a — + 0. 

Once the leading "IR divergent" cutoff effects 
have been canceled out, the use of the Symanzik 
expansion in the r.h.s. of eq. H53|l yields 



^2mq^ 

n„ Zj 



(54) 



where = |(^2|P±|7r±) l™"!^. In eq. JSl dots 
denote less dangerous "IR divergent" lattice arti- 
facts compared to those explicitely shown, as well 
as contributions of higher order in niq. In getting 
eq. (|54ll we have used the fact that the continuum 



limit of — (f2|S'gj^jj^|r2)^ at vanishing quark mass is 
E 7^ (see eq. 

From the above analysis it follows that, in the 
region where the series in eq. H54f) converges (i.e. 
at least down to masses where the order of mag- 
nitude inequality ~ niq > is still satisfied) , 
the squared mass of the charged lattice pion is 
linear in niq (up to very small deviations). Thus, 
vice- versa, we can imagine to use deviations from 
the established linear behaviour that are possi- 
bly seen at small niq as a workable criterion to 
determine the minimal value of niq at which sim- 
ulations can be performed before being set-off by 
discretization effects. 



5.4. Artifacts on hadronic energies and f-^ 

We wish to conclude by discussing some pe- 
culiar features concerning the magnitude of the 
O(a^) cutoff effects on hadron energies (in partic- 
ular masses) and on the pion decay constant. 
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5.4.1. Hadron energies 

In the language of the Symanzik expansion 
discretization artifacts on hadronic energies are 
described by a set of diagrams where at least 
one among the inserted / -Codd factors gets ab- 
sorbed in a multi-particle irreducible matrix ele- 
ment, with the consequence that it is not available 
for producing a pion pole. As a consequence, at 
fixed order in a, the most "IR divergent" lattice 
corrections to continuum hadronic energies con- 
tain one overall 1 / factor less than the leading 
"IR divergent" cutoff effects generically affecting 
correlators. For instance, at order d} the differ- 
ence between lattice and continuum energy of the 
hadron q;„ reads |6j (the label q specifying the 
three-momentum of the state |Q;„(q)) is omitted) 



AS„„(q) cx 

a" „ (f}|4k°(0))(7rO(0)a„|4|a„) 



(55) 



-Re- 



mt 



+ 0(a2) 



where O(a^) denotes "IR finite" corrections. It 
should be noted that Ai^Q^ (q)|(j2 is reduced to a 
plain O(a^) "IR finite" correction after anyone of 
the two "cures" described in Sect. 15.21 

Speciahzing (|55|) to the case of pions, one ob- 
tains the interesting result that the difference be- 
tween charged and neutral pion (square) masses 
is a finite O(a^) quantity even if the critical mass 
has not been set to its optimal value or the clover 
term has not been introduced. The reason is that 
the leading "IR divergent" contributions shown 
in are equal for all pions (as one can prove 
by standard soft pion theorems 0^1 )j hence cancel 
in the (square) mass difference. This conclusion 
is in agreement with detailed results from chiral 
perturbation theory (see refs. [20] and ^ODj as 
well as with the first numerical estimates of the 
pion square mass splitting in Mtm-LQCD |22|. 

5.4.2. Pion decay constant 

Data on a (quenched) computation of f-^, car- 
ried out in Mtm-LQCD by using the value of the 
critical mass obtained from the vanishing of the 
pion mass, show cutoff effects which tend to in- 
crease as the quark mass is lowered j^T]. In partic- 
ular data lie lower than the straight line extrapo- 
lation drawn from larger masses. This behaviour 



(called "bending phenomenon" in ref. TT) sets in 
at values of the bare quark mass around aAq^j-, , 
Aqcd ~ 200 MeV. A recent scaling test of Mtm- 
LQCD IS" furthermore shows that the "bend- 
ing phenomenon" is an 0{a?) cutoff effect with a 
magnitude which increases as m? ^ is lowered. 

In refs. I47I42I48] the lattice pion decay con- 
stant, was determined using the formula |i2i3'] 



2m, 



(56) 



from which is seen to be the ratio of two lattice 
quantities. It should be noted that the matrix 
element in the numerator of eq. (|56|l is affected 
by leading "IR divergent" lattice artifacts which 



are softened by an extra 



factor. This 



property can be traced back to the invariance of 
the correlator {P^{x)P^{y))\!^^ under axial ro- 
tations around the third iso-spin direction. One 
can, in fact, always think of having chosen the - 
0(a) - rotation angle such to bring the critical 
mass to its optimal value, at the price of shifting 
the bare quark mass by a term of order a^/niq. 
Also (m^±)^ is affected by leading "IR divergent" 
cutoff effects softened by an extra factor. The 
latter, however, drops out in the ratio {m^±)'^ /niq 
entering eq. fKHji . 

With anyone of the two cures described in 
sect. l5.2l lattice artifacts are reduced to terms that 
are only of the kind ^e>o ^^(a^/™?)^ 0- 

Indeed, a beautiful confirmation of the valid- 
ity of the analysis of "IR divergent" cutoff effects 
presented here comes from the fact that, when 
the critical mass is set at its optimal value, no 
bending effect is anymore visible in the data 
obtained in ref. |42l43j (see also P3|). 

6. Conclusions and outlook 

In this lecture we have reviewed the main prop- 
erties of Mtm-LQCD with the purpose of under- 
lining the features that makes it an appealing 
regularization for realistic simulations of QCD4 
at adequately small pion masses. Actually ex- 
plorative unquenched simulations of QCD4 with 
renormalized masses satisfying = md < riia < 
TOc have already started |48j . 

We have not discussed the very important is- 
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sue of meta-stabilities that on too coarse lattices 
appear to affect unquenched simulations and pre- 
vent reaching sufficiently light pions |21l48j . at 
least in the Singleton-Sharpe scenario [19j. 

In the language of lattice xPT the strength 
of meta-stabilities and the magnitude of the at- 
tainable minimal pion mass are features that are 
both controlled by the magnitude of the coeffi- 
cient (called C2 in refs. |19I2U| ') which in the chi- 
ral Lagrangian multiplies the term (Tr[S -I- S^])^ 
describing O(a^) lattice artifacts in correlators. 
Ideally, one thus would like to have C2 = as 
a function of g^. This situation is not peculiar 
of Mtm-LQCD. Similar features arise for generic 
values of the twisting angle and even for plain 
(clover) Wilson fermions a non- vanishing C2 might 
turn out to be problematic for simulations at fixed 
a if too small pion masses are taken |19j . 

Comparing xPT analysis with results from 
the Symanzik LEE A approach, one can prove 
that C2 is proportional to the matrix element 
a'^{Tr°{0)\e^-'°'\Tr°{0)), where e^~^' is the chirally 
breaking piece of (see eq. (03)- This observa- 
tion suggests that a possible way to enforce the 
condition C2 = is to modify the gauge action 
by adding to the standard plaquette term fur- 
ther dimension 6 operators. For instance, one can 
imagine adding &(.9o)^2xi, where 6(.g§) is an ad- 
justable coefficient and P2XI the minimal rectan- 
gular plaquette. This term, through the mixing 
induced by the chirally breaking twisted Wilson 
term, will modify, among others, all dimension 6 
operators of the Symanzik LEEA, including those 
that break chiral invariance. Since C2 also controls 
the square mass difference between charged and 
neutral pion |2UllUlllj . Am^ = rn^± — rn^o oc C2, 
the idea to ease Mtm-LQCD simulations is then 
to fix 6((7o) as a function of by imposing the 
condition Am J — 0. 
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